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Abstract. We prove that the expectation value of the index func- 
tion if(x) over a probability space of injective function / on any 
finite simple graph G = (V, E) is equal to the curvature K(x) 
at the vertex x. This result complements and links Gauss-Bonnet 

J2xev K ( x ) = X(G) and Poincare-Hopf J2 x ev = X(Q which 
both hold for arbitrary finite simple graphs. 



1. Introduction 

For a general finite simple graph G = (V,E), the curvature at a 
vertex x is defined as the finite sum 



k=0 

where Vk(x) is the number of K^+\ subgraphs in the sphere S(x) at 
a vertex x and V-\{x) = 1. With this curvature, the Gauss-Bonnet 
theorem [6] 

holds, where x{G) = YlT=o(~^-) hv k is the Euler characteristic of 
the graph, and where is the number of K^+i subgraphs of G. For 
example, if G contains no tetrahedral subgraph JQ, then each sphere 
S(x) lacks triangular subgraphs and K(x) = 1 — Vo(x)/2 + l / i(x)/3 and 
x(G) = |^| — \E\ + |T|, where T is the set of triangular subgraphs K 3 
of G. For an injective function / on the vertex set V, the index at a 
vertex is defined as the integer 

v(X) = 1 - x{S~(x)) , 
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where S~(x) = {y G S(x) \ f(y) < f(x) } is the exit set of the unit 
sphere S(x) with respect to the gradient field of / and where x{H) is 
the Euler characteristic of a subgraph H of G. The index if{x) is a 
discrete version of the Brouwer index for gradient vector fields and 
satisfies the discrete Poincare-Hopf theorem [7] 

53 1/0*0 = x(G) ' 

a result which holds for arbitrary simple graphs. Poincare-Hopf gives 
a fast way to compute the Euler characteristic because the subgraphs 
S~(x) are in general small. This allows to compute x{G) f° r ran_ 
dom graphs with hundreds of vertices, where counting cliques would 
be hopeless. 

Since Poincare-Hopf works for all injective /, also the symmet- 
ric index j f (x) = [i f (x) + i_/(z)]/2 satisfies J2xevJf( x ) = x(G). 
For cyclic graphs G, jf{x) is zero everywhere and agrees with cur- 
vature K(x). For trees, the symmetric index satisfies jf(x) — 1 — 
deg(x)/2 = 1 — |Vo(x)|/2 which adds up to 1 for connected trees. For 
graphs in which every unit sphere is a cyclic graph like the icosahe- 
dron, we have x(S + ( x )) = x(S~( x )) an d jf — if- The curvature 
K(x) is then 1 - \V Q (x)\/2 + |Vi(x)|/3 = 1 - \3(x)\/6 and the in- 
dex is if(x) = 1 — (^"(x)! = 1 — Sf(x)/2, where Sf(x) is the number 
of sign changes of / on the cyclic graph S(x). A small computation 
shows that in that particular case, the integral over all Morse functions 
E[s(x)] = |S'(x)|/3 and that E[l — s(x)/2] = 1 — \S(x)\/6 agrees again 
with curvature. This special case of the index expectation result led us 
to the more general result proven here. 

To keep this paper self contained, the proofs of the Gauss-Bonnet 
and Poincare-Hopf results are attached in an appendix. These gen- 
eral results become more geometric when dealing with graphs which 
are triangularizations of manifolds. In that case, Gauss-Bonnet is a 
discretization of Gauss-Bonnet-Chern and Poincare-Hopf is a discreti- 
sation of the analogue classical result in the case of gradient fields. 
In the continuum, for Riemannian manifolds, Euler curvature is only 
defined for even dimensional manifolds. This paper is a step towards 
proving that for odd dimensional graphs the curvature is always zero, 
something we know only in dimensions 1 and 3 so far. In an upcom- 
ing paper, using further developed techniques initiated here but using 
geometric assumptions on graphs like that unit spheres share proper- 
ties of the continuum unit spheres in d dimensions, we will prove that 
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for odd dimensional geometric graphs, the symmetric index jf(x) is 
zero everywhere. This matches the continuum case, where for Morse 
functions / the Brouwer index at a critical point is if(x) = (— l) m w 
where mix) is the Morse index, the number of negative eigenvalues of 
the Hessian matrix H(x) at the critical point x. In odd dimensions, 
of course = [if(x) + i_/(x)]/2 = at every critical point imply- 

ing immediately Poincare's result that odd- dimensional manifolds have 
zero Euler characteristic. We still are in search for continuum analogue 
of Theorem Q. The technical difficulty is to find a natural probability 
space of C 2 Morse functions on a compact Riemannian manifold. This 
is not a problem in the case of graph as we will see in the next section. 

2. Index expectation 

We first define the probability space of injective functions on the 
vertex set V of the graph G. Denote by n the order of the graph, the 
number of vertices in V . 

Definition. Let Q C [—1, l] n be the subset of all injective functions 
on V taking values in [—1, 1]. This is a n-dimensional Lebesgue space. 
We assume that Q is equipped with the product Lebesgue measure P. 
This means that P[{/ | f(x) G [a,b]}) = (b - a)/2 if -1 < a < b < 1 
and that the random variables X v (f) = f(v) giving the function values 
on the vertices v E V are independent and identically distributed. The 
injective functions are the complement of a union E of hyper surfaces 
in [—1, l] n and have full measure. Denote by E[i/(x)] the expectation 
of the index if(x) at the vertex x G V of / G f2 in this probability 
space (Q, P). 

In order to prove the main theorem, we need an excursion to percola- 
tion theory (We do not look at classical problems but for background, 
see O [3]), in particular site percolation, where vertices of a graph S 
are killed with a certain probability: given a background graph S and 
a fixed graph H, denote by vh the number of times the graph appears 
embedded in S. Now switch off vertices and edges connecting them 
in S independently from each other with probability p. Call v p H the 
expected number of graphs which appear now. It depends on S and p 
as well as H. We will see however that v p H dp only depends on H . In 
our case, we need the situation when H is the /c-dimensional simplex, 
the complete graph with k + 1 vertices. 

Denote by Vk(x) the number of H = Kk+i subgraphs in the sphere 
S(x) and by the number of K k+ i subgraphs in the exit set S~(x), 
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the subgraph of S(x) generated by vertices y where f(y) < f(x). Let vh 
denote the number of simplices K^+i which appear as subgraphs in S. 
We can look at the Erdos-Renyi probability space [1] of all subgraphs 
of S, where each vertex is included with probability p and the subgraph 
is the graph generated by these vertices. Let v p H the expected number 
of ^-dimensional simplices in the decimated subgraph of S. Obviously 
v v H < v ki but how much? Computing the expectation v p k of the survival 
rate depends on S and p But if p is chosen randomly too at first and 
each vertex is deleted with probability p, the survival rate only depends 
on the order of the clique and not on the graph: 

Proposition 1 (Clique survival for site percolation). 



Proof. The result is true if all the K^+i graphs in S are disjoint be- 
cause the survival of a single isolated simplex with k + 1 vertices is 



To prove the result in general, we decorrelate the situation by split- 
ting vertices: pick a vertex v where at least two such Kj. + i subgraphs 
Hi,H2 intersect. Replace v with 2 verticeswi , v 2 and place the edges 
to Hi with v\ and edges to H 2 with v 2 . Distribute the other edges 
originally intersecting with v arbitrarily with V\ or v 2 . To show that 
Jq 1 E p [vf] dp does not change when passing to the larger probability 
space, we compare the case before and after splitting: before splitting, 
v appears with probability p and contributes to Hi and H 2 . This gives 
Jq 2p 2 dp = 1. After splitting, both v.- t appear together with prob- 
ability p 2 , exactly one appears with probability 2p(l — p) and none 
appears with probability (1 — p) 2 . This also leads to a contribution 
2 JqP 2 dp + 1 Jq 2p(l —p) dp + Jq(1 —p) 2 dp = 1. We repeat like this 
with other intersection points of Hi and if 2- After all the correlations 
between Hi,H 2 are unlocked, we have a situation where the two sim- 
plices are independent and where the expectation value is the same as 
before. Now proceed with any other pair of simplices K^+i in the same 
way. □ 

The result means for k = that half of the points survive and for 
k = 1 that 1/3 of all edges are expected to survive. We ran Monte 
Carlo simulations with random host graphs G which is fixed over the 
experiment, where each vertex is knocked off with probability p. Ap- 
plying m such disaster experiments, each time starting fresh with the 
same G and then repeating the experiments for various p and averag- 
ing over disaster severeness p confirms the result remarkably well with 




J p k+1 dp= 1/0 + 2). 
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errors of the order 1/m and smaller simplices Kk+i like k = 2,3, 4. 

This result is actually more general. The clique graph K^+i can be an 
arbitrary pattern graph H . The result holds both for site and bond 
percolation situations. For site percolation catastrophes, the nodes 
are killed with probability p, in bond percolation catastrophes, the 

edges are broken with probability p. For an arbitrary background host 
graph 5* and any fixed pattern graph H, the expected decimation 
rate for the number of patterns H occurring in S is l/(ord(iJ) + 1) for 
site disasters and l/(size(iJ) + 1) for bond disasters. These network 
stability results are remarkably universal: they are independent of the 
background graph S. They can serve as "rules of thumb" if one has no 
a priory idea about the disaster strength p. 

Corollary 2 (Averaging equation). For every vertex x G V and all 
k>0, 

Proof. Look at a central vertex x connected to other vertices z%. We 
want to apply the previous lemma for S = S(x). Because / takes 
values in [—1, 1], we can assume f(x) = — l + 2p with p G [0, 1]. Having 
f(x) fixed like that, we get a random site percolation problem in the 
sphere S(x), where each vertex y G S(x) appears with probability 
p independently of each other. The expected number VjT(x) of k- 
dimensional simplices Kk+i divided by the number Vk(x) of simplices 
in S(x) is by the previous lemma equal to l/(k + 2) after we integrate 
over p. □ 

For k = 0, we have E[Iq + (x)] = E[V^ _ (x)] because the probability 
space is invariant under the involution / — > —f so that E|V ~(a:)] = 
Vq(x)/2 follows. For k — 1, the averaging equations are E[V r 1 ~(x)] = 
Vi(x)/3 and E[V{~ (x)} = E[y i + (a;)] = E[Wi(x)], where W x {x) is the set 
of vertices connecting vertices from V-^(x) to V{~(x). 

Here is the main result: 

Theorem 3 (Index expectation is curvature). For every vertex x, the 
expectation ofif(x) is K(x): 

E[i f (x)} = K{x) . 
Proof. We the averaging equation where k is replaced by k — 1 
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to see 



E[l- X (S-(x))} 



l-^(-l) fc E[^W] 

k=0 

oo 



fc=0 

oo 



(k + 2) 



fc=l 



fc vfc-l(s) 

(k+1) 



fc=0 



; T4-i(x) 
(* + l) 



□ 



Remark. This gives a new proof of the discrete Gauss-Bonnet result 
from 

$>(*) = x(G) 

simply by taking expectation. But unlike in the continuum, where 
Gauss-Bonnet-Chern is more difficult to prove (see e.g J2JH|), the dis- 
crete Gauss-Bonnet is easy to prove directly. As demonstrated in the 
Appendix, Gauss-Bonnet for graphs is even more direct than Poincare- 
Hopf. We expect that for compact Riemannian manifolds, a new prob- 
abilistic link between Poincare-Hopf and Gauss-Bonnet will allow to 
simplify the proof of the later considerably in higher dimensions. In 
the continuum, Poincare-Hopf is orders of magnitudes less complex 
than Gauss-Bonnet-Chern because it is part of differential topol- 
ogy, not needing any Riemannian metric while Gauss-Bonnet is part 
of differential geometry which uses more structure on the manifold 
M. It is the probability space on Morse function which will add part 
of the Riemannian structure on M, enough to get curvature. In the 
continuum, there are various probability spaces which are good candi- 
dates to represent curvature as index expectation. They all appear to 
work for compact two-dimensional surfaces. 
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Appendix 



Here are the proofs of Gauss-Bonnet [6] and Poincare-Hopf [7] for 
simple graphs G = (V, E) with consolidated notation. For Mathe- 
matica code, see [91 [10]. More Mathematica code illustrating all the 
probabilistic aspects proven in this paper and [8] will become demon- 
strations too. 

The first lemma generalizes Euler's handshaking lemma Ylxev ^o( x ) = 



Lemma 4 (Transfer equations). Ylxev Vk-ii^) — (k + l)v^. 

Proof. We can interpret I4_i(x) as the ^-degree of a vertex v, the 
number of fc-simplices Kk+\ which contain v. The sum over all k- 
degrees ^2 k deg k (x) is k + 1 times the number v k of /c-simplices K k +i 



Here is an other more pictorial proof of the transfer equations: draw 
and count handshakes from every vertex to every center of any k- 
simplex in two different ways. A first count sums up all connections 
leading to a given vertex, summing then over all vertices leading to 
J2xev ^4-i A second count is obtained from the fact that every 
simplex has k + 1 hands reaching out and then sum over the simplices 
gives {k + l)vk handshakes. 

Theorem 5 (Gauss-Bonnet). ^] y K(x) =x(G). 
Proof. By definition of curvature, we have 



Since the sums are finite, we can change the order of summation. Using 
the transfer equations (jlj), we get 



Given an injective function / on V and a vertex x G V, we can 
look at the set Wk(x) of all k simplices in the sphere S(x) for which at 
least one vertex y satisfies f(y) < f(x) and an other vertex z satisfies 



2vi. 



in G. 



□ 





□ 
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Proof. For each of the v k +i simplices K k+ 2 in G, there are k vertices x 
which have neighbors in K k+2 with both larger and smaller values. For 
each of these k vertices x, we can look at the unit sphere S(x) of v. 
The simplex K k+2 defines a k- dimensional simplex K k+1 in that unit 
sphere. Each of them adds to the sum Y^ x ^v which consequently 

is equal to kv k +i- □ 

Lemma 7 (Index stability). The index sum &v if(x) is independent 
off- 

Proof. The proof is a deformation argument. Fix a vertex x and change 
the value of the function f(x) such that a single neighboring point 
y G S(x), the value f(y) — f(x) changes sign during the deformation. 
Without loss of generality we can assume that the value of / at x has 
been positive initially and gets negative. Now S~(x) has gained a point 
y and S~(y) has lost a point. 

To see that x(S~(x)) + x(S~(y)) stays constant, we check this each 
individual simplex level and show V k (x) + V k {t) stays constant, where 
V k (x) denotes the number K k+ i subgraphs of S(x) which connect 
points within S(x) ± . Since i(x) = 1 — J2 k (— ^-Y^k i. x ) ■> ^ ne lemma 
is proven if V^~(x) + V^~(x) stays constant under the deformation. Let 
Uk(x) denote the number of Kj. + i subgraphs of S(x) which contain y. 
Similarly, let Uk(y) the number of K k+1 subgraphs of S(y) which do 
not contain x but are subgraphs of S~(y) with x. The sum of K k+1 
graphs of S~(x) changes by U k {y) — U k {x). When summing this over 
all vertex pairs x, y, we get zero. □ 

Theorem 8 (Poincare-Hopf). XLev = x(G) 

Proof. The number of /c-simplices V k ~ (x) in the exit set S~ (x) and the 
number of fc-simplices V k + (x) in the entrance set S + (x) are comple- 
mented within S(x) by the number W k (x) of k simplices which contain 
both vertices from S~(x) and S + (x). By definition, V k (x) = W k (x) + 
V k + (x) + V k ~(x). By the index stability lemma ([7]), the index if(x) is the 
same for all injective functions / : V — > R. Let x'(G) = ^ iel/ if(i). 
Because replacing / and — / switches S + with S~ and the sum is the 
same, we can prove 2v — J2xev x(S + (x)) + x(S~(x)) = 2x'{G) in- 
stead. The transfer equations Lemma (j4j) and intermediate equations 
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Lemma (jBD give 

oo 

2 X '(G) = 2^ + ^(-l) fe ^(V;-(x) + ^ + (x)) 



[1] 

[2] 

[3] 
[4] 

[5] 
[6] 

[7] 

[8] 

[9] 



fc=o xev 

DC 



k=0 xGV 

oo 

2^o + ^(-l) fc [(A; + 2)v k+1 - kv k+1 ] 

k=0 

oo 

2t'o + ^(-l) fc 2^ = 2 X (G) . 



k=l 

□ 
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